Optical Multi-hysteresises and "Rogue Waves" in Nonlinear Plasma 
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An overdense plasma layer irradiated by an intense light can exhibit dramatic nonlinear-optical 
effects due to a relativistic mass-effect of free electrons: highly-multiple hysteresises of reflection and 
transition, and emergence of gigantic "rogue waves". Those are trapped quasi-soliton field spikes 
inside the layer, sustained by an incident radiation with a tiny fraction of their peak intensity once 
they have been excited by orders of magnitude larger pumping. The phenomenon persists even in 
the layers with "soft" boundaries, as well as in a semi-infinite plasma with low absorption. 

PACS numbers: 42.65.Pc, 42.65.-k, 42.65.Tg, 52.35.Sb 

Diverse wave-related phenomena exhibit a common 
critical behavior: a dramatic transition at a crossover 
point from a traveling wave in the underlying (semi- 
infinite) medium to a non-propagating, evanescent wave 
that carry no energy. The crossover occurs in optics at 
the angle of total internal reflection at a dielectric inter- 
face, at a laser frequency near either a plasma frequency, 
or a waveguide cut-off frequency, or band-gap edge of a 
material, including BEC; in quantum mechanics for elec- 
tron scattering at the energy close to a potential plateau, 
etc. It can be of great significance to nonlinear optics: a 
nonlinear refractive index can cause a phase-transition- 
like effect, since a small light-induced change may trans- 
late into a switch from full reflection to full transmission, 
resulting in a huge hysteresis. Predicted in [1] for non- 
linear interfaces, it was explored experimentally in [2], 
with an inconclusive outcome, with some of the experi- 
ments (including the latest [2c]) showing a clear hystere- 
sis, while [2b] showing none (see also below). 2D numer- 
ical simulations were not well suited then for modeling 
hysteresis; their very formulation excluded multivalued 
outcome by using single-valued boundary conditions. 

In this Letter we discover, however, that even the 
most basic, ID-case, reveals a large and apparently lit- 
tle known phenomenon of highly-multistable nonlinear 
EM- propagation and the emergence of trapped "rogue" 
(R) waves, with intensity exceeding the incident one by 
orders of magnitude. They may be released by plasma 
expansion, in a phase-like transition, as in a boiling liq- 
uid. Even a slight nonlinearity due to the most funda- 
mental mechanism - relativistic (RL) mass-effect of free 
electrons - suffices to initiate the effect. Multiple (up to 
hundreds) hysteretic jumps between almost full reflection 
and full transparency may occur as the laser intensity is 
swept up and down. We treat the problem here in the 
context of ideal plasma on account of recent interest in 
high-intensity laser-plasma interactions and fundamental 
nature of RL-nonlinearity, but all approaches are valid 
for other crossover problems. Temporal RL-solitons have 
been considered in detail in the literature [3]; close to 
those of underdense plasma are the so called Bragg or 
band-gap solitions [4] , including the ones in BEC [5] . The 



difference in this work is made by multi-hysteresises (and 
standing, unmoving quasi-solitons instead of propagating 
ones) due to self-induced retro-reflection. The remark- 
able new property is that for the same incident power an 
EM-wave can penetrate into a nonlinear material to dif- 
ferent depths - varying by orders of magnitude - depend- 
ing of the history of pumping. We assume a stationary, 
cio, or long pulse mode, and use only RL-nonlinearity in 
a cold plasma. While this model is greatly simplified vs 
various kinetic approaches, it allows us to keep the basic 
features necessary to elucidate new results, and have the 
theory applicable to other systems. 

Even a few-A-thick plasma layer can produce the ef- 
fect, so that the absorption of light (including nonlinear 
self-focusing) would not affect the propagation signifi- 
cantly. Furthermore, while abrupt boundaries contribute 
to the effect, they are not essential: a layer with "soft" 
boundaries still exhibits all the major features of the 
phenomenon. The phenomenon persists even for semi- 
infinite plasma with a small absorption, which also de- 
velops a strong retro-reflection. The " Sommerfeld condi- 
tion" (no wave comes from the " infinity" ) is to be revis- 
ited here: a wave is back-reflected deep inside the plasma 
and returns to the boundary. Imposed on a forward wave, 
it results in a semi-standing wave, trapped quasi-soliton 
R- waves, and multistability, same as in a finite layer. 
The energy accumulated in R- waves and excited free elec- 
trons, can then be released if plasma density reduces, as 
it may happen, e. g. in an astrophysical environment. 

The wave propagation here is governed by the same so 
called nonlinear Klein-Gordon-Fock equation: 

[v 2 - a 2 /(v9i) 2 ]^ = fcg/(M 2 )V; /(o) = i. (i) 

where /(|"0| 2 ) is a function responsible for nonlinearity. 
Here, a generic variable, if, could be a scalar (e.g. a 
wave function in RQM), or a field vector in EM-wave 
propagation, v is a scale velocity, (v = c for plasma and 
in RQM), and k^ 1 is a spatial scale of the problem, (for 
RQM, k — moc/h — 2it/Xc, where Ac is the Compton 
wavelength, and mo is the rest-mass of an electron); for 
plasma fco = w p //c, where co p i is a plasma frequency due 
to free-charge density, p e , and in the X-ray physics, Ui p h — 
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Eph/h is related to photo- ionization limit, E p h, of atoms. 
For an ^-monochromatic wave ip = £(r) e~ tuJt /2 + c.c, 
where £{r) is a complex amplitude, and using u 2 = \£\ , 
(1) is reduced to a nonlinear Helmholtz equation: 

V 2 £ + k 2 e(u 2 )£ = 0; e pl = 1 - uj 2 pl {u 2 ) / u: 2 (2) 

where k — uj/c, and cj 2 ; = 4irep e /m, m - electron mass, 
and £ = E/E n i with E n i being some characteristic non- 
linear scale; for a RL-mass-effect it is E r i = uimoc/e. 

In general, nonlinear e may have various origins: a 
varying ionization rate, plasma waves, ponderomotive 
force, etc. Assuming fully ionized gas, p e = const, and a 
circularly-polarized wave, £(£) (e x + ie y ) e~ lult /2 + c.c, 
that has very negligible high-harmonics generation and 
minimal longitudinal plasma waves excitation, the most 
basic remaining source of nonlinearity is a field-induced 
RL mass-effect of electron: m — mo7, with a relativis- 
ts factor 7 = \JT+ (p/m a c) 2 = \/l + u 2 [see e. g. [6]) 
where p is the momentum of electron, so that 

e r i — 1 - [f 2 7(w 2 )] _1 with v = uj / (ojpi)o (3) 

where {u) p i)o is a linear plasma frequency with m = mo. 
Since m = m{u 2 ), a single electron exhibits large hys- 
teretic cyclotron resonance predicted in [6a] and observed 
in [7]. The mass-effect has also became one of the ma- 
jor players of light-plasma interaction [3], e. g. in RL 
self-focusing, and in acceleration of electrons by the beat- 
wave and wake-field. The EM-propagation could also be 
accompanied by RL-intrinsic bistability [8]. 

In a ID-case, letting a plane EM- wave propagate in the 
z-axis, we have V 2 = d 2 /dz 2 . For a boundary between 
two dielectrics, a EM- wave incident from a dielectric with 
tin > under the angle 9 onto a material of tNL > 0, 
e in (2) is replaced by ei n [^NL(^) / 'fm(w)- sin 2 #]. For 
a mw waveguide with a critical frequency uj wg , e in (2) 
is replaced by e wg (l — aj^g/u 2 ). The crossover point is 
attained at e = 0. In this approximation, (2) reduces to 

£" + e((,u 2 )£ = 0, (4) 

where £ = kz, and "prime" denotes d/d(; in general, 
we do not assume e uniform in £-axis. In a weakly- 
nonlinear media one can break the field into counter- 
propagating traveling waves and find their amplitudes 
via boundary conditions. However, near a crossover 
point one in general cannot distinguish between those 
waves. To make no assumptions about the wave compo- 
sition, we represent the field using real variables u, and 
phase (eikonal), </>, as £ = u(()exp[i<p(Q]. Since £ is in 
general complex, while e = e(u 2 ), Eq. (4) is isomorphous 
to a 3-rd order eqn for u; yet, it is fully integrable in 
quadratures. Its first integral is a scaled momentum flux 
of EM-field, P = v?4> = inv. In a lossless media P is con- 
served over the entire space £ < oo, even if the medium is 
non-uniform, multi-layered, linear and/or nonlinear, etc. 



If a layer borders a dielectric of e = e ex at the exit, we 
have P = u 2 x ^Jt ex , where u 2 x is the exit wave intensity. 
Eq. (2) is reduced then to a 2-nd order equation for u: 

u" + u[e((,u 2 )-P 2 /u 4 } =0, (5) 

which makes an unusual yet greatly useful tool. By ad- 
dressing only a real amplitude, while using flux P as a 
parameter, (5) is nonlinear even for a linear propaga- 
tion, yet is still analytically solvable if a density p e is 
uniform across the layer {de/dC, = 0). A full-energy- like 
invariant of (5) is u' 2 /2 + U(u 2 ) = W = inv, with U = 

[jf e(u 2 ) d(u 2 )+ P 2 /u 2 }/2 , where u' 2 /2 is "kinetic", 
and U - "potential" energies. For a RL-nonlinearity (3), 
U(u 2 ) = (u 2 + P 2 /u 2 )/2 - Y){u 2 ) - l]/v 2 . Here W is 
a scaled free EM energy density of e-nonlinear medium 
[9] W = c[H 2 + J ed{E 2 )]/{2E 2 rl ), where H is magnetic 
field. If a layer exit wall is a dielectric, one has W — 
U(u 2 x ), since then u' ex = (see below). For a metallic 
mirror, W — u' ex j2, since now u ex — 0; and W = 
for an evanescent wave in a semi-infinite medium. The 
implicit solution for spatial dynamics of u in general case 
is found now as ( = J{2[W -[/(u 2 )]}" 1 / 2 ^ . 

Boundary conditions at the borders with linear di- 
electrics at the entrance, £ = 0, with t in , and at the 
exit, ( — d, with e ex , result in complex amplitudes of 
incident, £ in , and reflected, £ r fi, waves at £ = 0: 

£ in , rfl = [u± e; n 1/2 (P/u - m')]/2; (6) 

where "+" corresponds to £, n , and "— " - to £ r fu and 
the exit point, £ = d: u — £ ex = u ex ; u' — 0; and <f>' — 
y/e^. A condition P — corresponds to full reflection, 
resulting in either strictly standing wave, or nonlinear 
evanescent wave in a semi-infinite plasma, in particular 
in a " standing" soliton-like solution (see below) . Ife(u 2 = 
0) < 0, there arc no linear traveling waves; yet a purely 
traveling nonlinear wave may exist at sufficiently strong 
intensity u 2 = u 2 rv = const, such that e(u 2 rv ) > 0: 

£trv = u trv exp(i<//£), (j)' = ±^Je{u 2 trv ), (7) 

propagating either forward (+) or backward (— ). How- 
ever, if e(u = 0) < 0, it is strongly unstable. A non- 
periodic solution of (5) with P — is a nonlinear evanes- 
cent wave that forms a standing, trapped soliton. In 
low-RL case, one needs a small detuning from crossover 
point, 5 = v — 1 « 1, to attain the effect at low laser 
intensity, u 2 -C 1, so that the the dielectric constant (3) 
is Kerr-like and small: e T i rts —2S + u 2 /2, \e r i\ <1. A full 
solution of (5) with P = and u — > at £ — > oo yields 
then a standing soliton with a familiar intensity profile: 

u 2 = 85/ cosh 2 [{( - ( pk )\f26] (8) 

where the peak location ( p k is an integration constant. 
For an arbitrary frequency, v < 1, the soliton peak in- 
tensity is u 2 ol — 4(1 — v 2 )jv A , instead of 88 as in (8); 
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an unbound plasma, the solution of (5) with a traveling 
component, P > 0, is a spatially periodic and positively 
defined, with the intensity, u 2 ((), bouncing between two 
limits, ul x , and u 2 k . If P/16 « 5 2 « 1, we have 



u\ x < u 2 < u 2 pk 



8<5 + P/26 



(9) 
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FIG. 1: Hysteretic transmission of light through a plasma 
layer of thickness L. FT and NTR are full transparency and 
near total reflection limits. Points 1, 2, 3 mark a linear evanes- 
cent wave, 1-st, and 2-nd upper stable states respectively, and 
4 - a R-wave sustained by a very low pumping. Arrows in- 
dicate direction of jumps within the lowest hysteretic loop. 
Inset: spatial amplitude profiles of waves corresponding to 
points 1-4 in the main plot. 



i. e. the peaks are relatively large, u 2 
a train of well separated quasi-solitons nearly coinciding 
with a standing soliton (8) of the peak intensity u 2 k ~ 8(5. 
As P and u 2 n increase, they grow larger and closer to each 
other. The spatial period, A, of this structure is: 



A/A sa ln(166/V r P)/(2TrV26) 



(10) 



In a strongly RL case, P > 1, we have u 2 k /P w 
1 + (1 + P)~ 2 , and A w A/2, as for a standing, albeit in- 
hibited wave in free space, while traveling wave emerges 
dominant, resulting in self-induced transparency. 

Hysteretic jumps occur when either valley or peak of 
the intensity profile coincide with an entrance, £ = 0. 
The valley, uf t = u 2 xl marks an off-jump in Fig. 1, and 



the peak, Uq 



e(u 2 
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an on-jump. Suppose that in a layer 



pk 



) = (1 - v 2 )/{2 - v 2 ) > 0. When v 2 < 1/2, it is of L > \/(2irV26), the incident intensity u 2 n is ramped 



a strongly- RL soliton, u 2 ol 3> 1, and its peak narrows 
down to a half- wave: u 2 rts u 2 ol cos 2 (C — ( P k) at u 2 > 1. 

In a finite layer one has a mix of standing/evanescent 
and traveling waves, with v} min = u 2 x = P > 0. A full 
integration of (5) with nonlincarity e r ; « —2(5 + u 2 /2 
yields then elliptic integrals of imaginary argument of 
the first kind; more importantly, cq. (5) and its invari- 
ant avails themselves to detailed analysis. The numerical 
simulations are needed, however, to find a solution for 
(a) strongly- RL field [using (5) and its integrals], or (b) 
non-uniform plasma density in (5), or(c) plasma with ab- 
sorption [cq. (4)]. It then is found by an "inverse propa- 
gation" procedure, whereby we essentially back-track the 
propagation from a purely traveling exit wave back to the 
entrance. One sets first a certain magnitude of u 2 x = P, 
u' ex = at the exit, numerically computes an amplitude 
profile u(Q back to the entrance and incident and re- 
flected intensities u 2 x and u 2 ^ using (6), and then maps 
u 2 x and u 2 ^ vs incident intensity, u 2 n . A data set u 2 n (P) 
and v%fi(P) for any given P is found then with a sin- 
gle run, vs a so called multi-shooting commonly used in 
search of solution with conditions set at two boundaries. 
This provides a very fast numerical simulation vs multi- 
shooting; besides, the latter one is very unreliable when 
dealing with apriory unknown number of multi-solutions. 

For a fully-RL simulation with a, L ~ 10A, where L is 
the layer thickness, Fig. 1 show the emergence of large 
number, Nh s , of huge hysteretic loops of the transmission 
(same as in reflection, not shown here), which bounces 
between full transparency (near the points touching an 
FT line) to nearly full reflection (near the points touch- 
ing an envelope NFR). In general, Nh s — 0(L/X). In 



up from zero. When uf n < 28, Fig. 1. point 1, the 
amplitude is almost exponentially decaying, ( = 0, as 
u w 2ui n exp (—y/25£), i. e. is a nearly-linear evanescent 
wave, curve 1 in Fig. 1 inset; the layer is strongly refrac- 
tive, and the transmission is low. As u 2 n increases, the 
front end of that profile swells up, becoming a semi-bell- 
like curve, close to (8) with ( pk w 0. With further slight 
increase of pumping, it gets unsustainable, and the field 
configuration has to jump up to the next stable branch of 
excitation, whereby it forms a steady R-wave at the back 
of the layer. If after that uf n is pulled down adiabati- 
cally slow, the R-wave moves to the middle of the layer 
(Fig. 1. point 2, curve 2 in the inset). Finally, when it 
is exactly at the midlayer (Fig. 1. point 4, curve 4 in 
the inset), both valleys are at the borders of the layer, 
the pumping is nearly minimal to support an R-wave; 
below it the profile is unsustainable again, and the sys- 
tem jumps down to a regular nearly-evanescent wave and 
almost full reflection. 

At this remarkable point, the layer is fully transpar- 
ent, i. e. all the (very low) incident power is trans- 
mitted through, while a giant R-wave of peak intensity 
(ui n ) 2 k w 2(5 inside the layer is sustained by a tiny in- 
cident power, (uirOmin- ^ L^/26 / N \ > 1, the contrast 
ratio - essentially a nonlinear resonator's finesse, Q, is 



Q = 



(u 
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•n) p k 
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> 1 



(11) 



where N is the number of R- waves in a layer; the one with 
N = 1 occurs after the first jump-up. In the example for 
Fig. 1 (6 = 0.02, L = 10A), Q - 10 7 . In semi-infinite 
plasma, Q is limited by absorption, see below. It also 
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FIG. 2: Hysteretic absorption of light in a semi-infinite 
plasma layer with absorption a. FA is a full absorption limit. 
Points 1, 2, 3 mark a linear evanescent wave, f-st, and 2-nd 
upper stable states respectively. Arrows - the same as in Fig. 
f. Inset: intensity profile; points and verticals 1-3 indicate 
locations of the plasma boundary for the respective points in 
the main plot; <'s show direction into plasma layer. 



decreases as N increases; the field profile for TV = 2, is 
depicted in Fig. 1, point 3, curve 3 in the inset. 

Only half of multi-steady-states are stable; the stability 
condition is that the EM-energy density increases with 
the pumping, i. e. dW/d(uf n ) > 0, which also coincides 
with the condition d(u 2 x ) j d(iq n ) > 0, similar to [1]. 

One can view R-waves at a N-th stable branch as a 
N-th order mode of a self-induced resonator, with full 
transparency points marking the resonance. The mirror- 
like sharp boundaries of plasma layer enhance the reso- 
nances [10], but do not constitute necessary condition. 
Our numerical simulations using (4) showed that a layer 
with "soft" shoulders making ~ 50% of the entire layer 
length, still exhibits a few hysteresises, and a large num- 
ber of self-induced resonances. 

The fundamental manifestation of the phenomenon 
transpires in a semi-infinite plasma. Only two kind of 
waves [1] in a lossless case satisfy then the Sommerfeld 
condition - no wave "comes back" from £ — > oo - a 
traveling (7), du/dQ — > 0, and an evanescent (8) wave, 
u — > 0. Our investigation of (2), to be published else- 
where, showed that the wave (7) is unstable both in 
2D&3D-propagation - and, of (1) - in temporal domain 
in ID-case. However, using (4), one can show that even 
a steady ID-wave (7) is unstable against small absorp- 
tion by replacing v 2 in e p i (3) with v 2 {l + ia), where 
a = (cot) is an absorption factor, and r is an electron 
momentum relaxation time (typically a -C 1). A condi- 
tion for a hysteresis to emerge is then a < a cr rts eS [11], if 
S <C 1. Near a ~ a cr , a jump- up occurs at uf n w nS. The 
ratio (11) to sustain a single R-wave is limited now by 
Q rts S/a, and can still be huge. Hysteresises in reflection 
at a = 10~ 3 are shown in Fig. 2, and the intensity profile 
for u 2 n > 25 - in the inset. An initially traveling wave 



develops oscillations due to rising standing wave, which 
eventually becomes a train of trapped R-waves, the last 
one being a quasi-soliton close to (8), and then vanishes 
exponentially. Reducing a pushes that last R-wave fur- 
ther back, but does not extinguish retro-reflection from 
the R-wave train at the crossover area deep inside plasma, 
keeping the condition u — > at ( — »■ oo. 

Lab observation of the phenomena in plasma could be 
set up with e. g. jets of gas irradiated by a powerful C0 2 
laser, with a gas density controlled to reach a crossover 
point at Xco 2 ~ 10/um; the recent experiments [12] may 
actually indicate R-waves released by plasma expansion. 
This process may be naturally occurring in astrophysical 
environment in plasma sheets expelled from a star (e. g. 
the Sun); part of the star's radiation spectrum below the 
initial plasma frequency is powerful enough to penetrate 
into the layer and be trapped as R-waves. When the layer 
expands, they can be released as a burst of radiation. It 
is also conceivable that the R-waves trapping and conse- 
quent release may be part of the physics of ball-lighting 
subjected to a powerful radiation emitted by the main 
lighting discharge in mw and far infrared domains. The 
R-waves might be used e. g. for laser fusion to deposit 
laser power much deeper into the fusion pallets; or for 
heating the ionosphere layers by a powerful rf radiation. 

In conclusion, optical multi-hysteresises may emerge 
in a plasma near critical plasma frequency due to fun- 
damental relativistic mass-effect of electrons. They may 
result in huge trapped, or standing rogue waves with the 
intensity greatly exceeding that of pumping radiation. 

This work is supported by the US AFOSR. 



[1] A. E. Kaplan, JETP Lett. 24, 114 (1976); also Sov. 

Physics JETP 45, 896 (1977). 
[2] P. W. Smith, et al, Appl. Phys. Lett. 35 846 (1979); 

also IEEE JQE, 17 846 (1981); P. W. Smith and W. J. 

Tomlinson, IEEE JQE, 20 30 (1984); S. De Nicola, et al, 

Appl. Phys. B 49, 441 (1989). 
[3] G. A. Mourou, T. Tajima, S. V. Bulanov, Rev. Mod. 

Phys., 78, 309 (2006) 
[4] W. Chen and D. L. Mills, Phys. Rev. Lett. 58, 160 (1987), 

J. E. Sipe and H. G. Winful, Opt. Lett. 13, 132 (1988), 

D. N. Christodoulides and R. I. Joseph, Phys. Rev. Lett. 

62, 1746 (1989); B. J. Eggleton, et al, ibid, 76, 1627 

(1996); D. Mandelik, et al, ibid, 90, 053902 (2003). 
[5] O. Zobay, S. Potting, P. Meystre, E. M. Wright, Phys. 

Rev. A 59, 643 (1999) 
[6] A. E. Kaplan, Phys. Rev. Lett. 48, 138 (1982); ibid, 56, 

456 (1986); A. E. Kaplan and P. L. Shkolnikov, Phys. 

Rev. Lett. 88, 074801, (2002). 
[7] G. Gabrielse, H. Dehmelt, and W. Kells, Phys. Rev. Lett. 

54, 537 (1985). 
[8] R. J. Noble, Phys. Rev. A 32, 460 (1985); B. M. Ashki- 

nadze and V. I. Yudson, Phys. Rev. Lett. 83, 812 (1999); 

G. Shvets, ibid, 93, 195004 (2004). 
[9] B. Y. Zeldovich, Brief Comm, Lebedevlnst, 2, 20 (1970). 



■5 



[10] Weakly nonlinear Fabri-Perot resonator has been used to 
first demonstrate optical bistability (first-order hystere- 
sis), H. M. Gibbs, S. L. McCall, and T. N. C. Venkatesan, 
Phys. Rev. Lett. 36, 1135 (1976); H. M. Gibbs, Optical 
Bistability (Academic Press, New York, 1985). 

[11] This suggests an explanation of the hysteresis lacking in 



[2b] : an artificial nonlinearity used there was due to small 
dielectric particles suspended in a liquid and had large 
dissipation because of strong light scattering, vs the ex- 
periments [2a] and [2c] with nearly transparent fluids. 
[12] P. L. Shkolnikov, private communication. 



